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Abstract 

We construct a model of quintessence in string theory based on the idea of axion 
monodromy as discussed by McAllister, Silverstein and Westphal [1]. In the model, 
the quintessence field is an axion whose shift symmetry is broken by the presence of 
5-branes which are placed in highly warped throats. This gives rise to a potential for 
the axion field which is slowly varying, even after incorporating the effects of moduli 
stabilization and supersymmetry breaking. We find that the resulting time depen- 
dence in the equation of state of Dark Energy is potentially detectable, depending 
on the initial conditions. The model has many very light extra particles which live 
in the highly warped throats, but these are hard to detect. A signal in the rotation 
of the CMB polarization can also possibly arise. 
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1 Introduction 



Observational evidence shows that the universe is accelerating and therefore must be dominated, 
at cosmological scales, by a mysterious form of energy called Dark Energy. Understanding its 
nature is a central challenge that faces us today. The leading candidate for dark energy is the 
cosmological constant. It is consistent with observations. Theoretically, we now know that a 
positive cosmological constant can arise from a consistent theory of quantum gravity like string 
theory and this gives us greater confidence in the idea being correct. 

Another possibility for dark energy is quintessence, see [2H6] for reviews. In this case the 
vacuum energy is not a constant but instead slowly relaxes due to the evolution of a scalar field. 
This makes the resulting equation of state for dark energy time dependent. Observations in the 
coming decade will attempt to determine the equation of state of dark energy and, hopefully 
will allow us to decide between these possibilities. 

In this paper we will construct a model of quintessence in string theory. One can view this 
as a sort of "theoretical test" of this idea. If the idea fits in with string theory one would have 
greater confidence in it, just as for the cosmological constant. One can also hope that such 
constructions might ultimately lead to some interesting constraints which can then be tested 
observationally. 

The idea behind quintessence is quite similar to that in inflation. However, an important dif- 
ference, which makes the construction of models of quintessence considerably more challenging, 
is that the energy scale for quintessence is of order E 



A~l(r 3 eV (1.1) 

and therefore is much smaller than the scale of supersymmetry breaking. This leads to two 
important issues which any model of quintessence must confront. The first is to ensure that the 
potential for the quintessence field has a scale of order A and not of order the SUSY breaking 
scale; this issue is related to the cosmological constant problem. The second is to have a potential 
for the quintessence field which meets the slow-roll conditions despite the relatively high scale 
of supersymmetry breaking. In particular this requires that the mass of the quintessence field 
is of order the Hubble constant today, H ~ 10 -33 eV, and therefore is much smaller than the 
SUSY breaking scale, which must be at least M$b ~ 1 — lOTeV and typically is much higher 
§. Now, it is a famous fact, responsible for the hierarchy problem in the standard model, that 
scalar fields tend to have their mass driven up to the SUSY breaking scale. Ensuring that the 
quintessence field has a mass m ~ H so that the ratio itl/Msb is at least 10 -45 if not smaller is 



: We will use A throughout in the following discussion to denote this energy scale. In particular A will not 
denote the cosmological constant in this paper. We hope this does not lead to any confusion. 

2 Here the supersymmetry breaking scale refers to the underlying scale at which the symmetry breaks, related 
to the gravitino mass by, m g ~ Mg g /Mpi. 
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therefore not easy. 

In this paper we will, to begin, ignore the first question, related to the cosmological constant 
question, and focus on the construction of a model which meets the second challenge of ensuring 
that the slow-roll conditions are met despite the high scale of SUSY breaking. We will return to 
a discussion of the cosmological constant question in §6. Within the context of a model of the 
type we construct, we will argue that anthropic considerations could well result in the energy 
density in dark energy being shared between the cosmological constant and the quintessence 
field in a reasonably equitable fashion. As a result an evolving equation of state might well be 
more natural from the point of view of anthropic considerations in such a model. 

Before proceeding it is useful to phrase the requirement for a small mass for the quintessence 
field in terms of restrictions on higher dimensional operators in the low-energy effective field 
theory. Taking M$b ~ 1 TeV, and the mass m ~ H, one finds that for the current value of 
H ~ l(T 33 eV, 



This shows that upto dimension 10 operators which could contribute to the mass have to be 
suppressed to construct a satisfactory model. For a higher scale of SUSY breaking the sup- 
pression would have to be even more severe. The requirements of slow-roll are therefore indeed 
quite restrictive §. This estimate also illustrates why the construction of a quintessence model 
is sensitive to Planck scale physics and can therefore be sensibly carried out only within a UV 
complete theory of gravity such as string theory. In fact in the model we construct the canoni- 
cally normalized quintessence field typically undergoes an excursion of order the Planck scale in 
the course of its evolution and the potential needs to be flat for this whole range in field space. 
Needless to say, ensuring this while working only in an effective field theory would be extremely 
difficult. 

Qualitatively, there are two possibilities within string theory for a quintessence field [7]. It 
could be a modulus like the dilaton or the overall volume which runs off to infinity in field 
space. In such a situation the coupling constants of the standard model, like the fine structure 
constant or the four dim. gravitational constant, would also vary with time and one would 
also need to ensure that this variation is small enough. Such models would be observationally 
very interesting but are quite challenging to construct. Instead, we will opt here for a second 
possibility which seems easier to implement. In our case the field will be an axion, which is a sort 
of Goldstone boson, that runs to a finite point in field space where its potential is minimized. 
We will take this finite point to be at the origin. 

3 The sensitivity to UV physics would be significantly reduced if we were sure that a global symmetry was left 
intact by it and only broken at low energies giving rise to the mass. However, in general Planck scale physics 
is expected to break global symmetries, so to be sure that the global symmetry is approximately unbroken one 
needs the full UV theory. We thank L. McAllister for related discussion. 
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It is well known that axion fields often arise in string compactifications. In our construction, 
which is in IIB string theory, the axion of interest will arise from the RR sector. A general 
argument then ensures that there is a shift symmetry, under which the axion shifts by a constant, 
which is unbroken to all orders in the quantum loop expansion. This feature will be very helpful 
in ensuring that the resulting potential is varying slowly enough and can be isolated from the 
effects of supersymmetry breaking. 

The shift symmetry will have to be broken of course, at some level, to generate a potential 
for the axion. A natural possibility to consider is that the breaking arises from non-perturbative 
effects like instantons. It turns out though that this does not lead to an acceptable model for 
quintessence @. In this case the resulting potential is a periodic function of the axion. The 
slow-roll conditions then require the axion to have a decay constant of order the Planck scale or 
higher whereas requiring the energy density in the axion to be of the right magnitude gives an 
axion decay constant about two orders of magnitude smaller than the Planck scale [HIIS] - 

One therefore needs to explore other ways to break the shift symmetry. It turns out that 
the shift symmetry can also be broken in the presence of branes §. By placing these branes in 
highly warped regions of the compactification one can make this breaking small §. This idea 
is referred to as "axion monodromy", and was developed in [IJQIBCEI] to construct a model of 
chaotic inflation. The resulting potential one gets in this case is no longer periodic, in fact it 
is approximately linear in the axion . As a result, we will find that the slow- roll conditions 
and other requirements for quintessence can be met even with an axion whose decay constant 
is smaller than the Planck scale. 

Following [1] (see also [H]), we will add a pair of NS5-brane and NS5-brane in highly warped 
regions of the compactification to break the axion shift symmetry. The axion, which is the zero 
mode of the RR two-form C2 that arises due to a non-trivial two-cycle S, induces a D3-brane 
charge on the 5-brane which wraps the two-cycle E. It also induces D3-brane charge on the 
anti 5-brane 0. This results in additional 3-brane tension being induced on both the 5 and anti 
5-branes which depends on the axion. The resulting potential turns out to be linear in the axion, 
for large enough values of this field. 

To examine the effects of SUSY breaking and moduli stabilization on the axion potential we 

4 With 0(100) axions rolling coherently one can get a workable model [8]; at the very least this is an invitation 
to seek other possible solution. 

For example, a Dp brane in the presence of the NS two-form B 2 acquires an induced D(p — 2) form charge 
and additional tension. This breaks the B2 shift symmetry although a symmetry still remains under which B2 
and the world volume gauge field F2 both shift keeping B2 — F% invariant. Similarly, the shift symmetry for C2 
is broken in the presence of NS5-branes. 

6 At first sight this way of breaking the shift symmetry might seem rather contrived. But such branes or related 
fluxes can in fact arise quite generically in flux compactifications. We thank E. Silverstein for a discussion on 
this point. 

7 The anti 5-brane is actually a 5-brane wrapping cycle £ but with opposite orientation. 
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will consider the model of KKLT [12J. We will find that an additional important contribution 
to the axion potential does arise from these effects. The axion affects the warp factor of the 
geometry and this in turn changes the volume of the internal space. Once the volume is stabilized 
this gives rise to an extra potential for the axion. The leading contribution of this type would 
in fact be unacceptably big, but this can cancel exactly between the 5-brane and the anti 5- 
brane which induce opposite 3-brane charge, provided the two throats where they are located 
are related by a Z2 symmetry. 

The subleading correction to the axion potential which survives cannot be calculated exactly, 
given our present understanding of flux compactifications, but an estimate which is adequate for 
the purposes of our discussion can be made following the discussion in [T3lTf5] . This estimate 
tells us that the subleading effect is still important enough to be the dominant contribution to 
the axion potential, but it too is suppressed by the warp factor at the bottom of the throats 
where the 5-branes are located §. By making this warp factor small enough one then finds that 
a satisfactory model of quintessence can indeed be constructed. 

The bottom-line then is that a workable model for quintessence in string theory, based on 
the idea of axion monodromy, can be constructed by placing a 5-brane and anti 5-brane pair 
in highly warped throats in the presence of an axion. When the canonically normalized axion 
field has a expectation value which is of order the Planck scale, the slow-roll conditions are met, 
and the quintessence field evolves with an equation of state which can deviate significantly from 
that for the cosmological constant. 

This paper is structured as follows. Some general considerations for models based on a non- 
perturbative potential and a linear potential are discussed in §2. The quintessence model we 
construct is then discussed in greater detail in §3. Additional terms which can arise once the 
dynamics of moduli stabilization is included are discussed in §4. §5 contains a more complete 
look at the final model including all important terms in the potential, and discusses the resulting 
energy scales in more detail. Some general features of the model are discussed in §6. These 
include the many very light particles which arise due to the highly warped throats in which the 
5-branes are placed, the absence of tracker behavior and its relation to the cosmological constant 
issue, and a potentially interesting signal due to the rotation of the E mode of the CMB into 
the B mode. We end with conclusions in §7. 

Let us end this introduction by commenting on some additional papers of relevance. A model 
of quintessence based on a linear potential was considered in [16] , and was explored further 
in [T71[l~8] [§. The idea that quintessence might arise from many axion due to the effects of flux 
induced mixing was explored in |19j . The possibility of many light axions and the constraints 

8 Although the suppression is by two powers of the warp factor, whereas the 3-brane tension contribution 
mentioned above is suppressed by four powers of this factor. 

9 We thank A. Linde for bringing these references to our attention. 
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that can be placed on them was discussed in [201 121] ■ Finally [22J, which appeared recently, 
discusses how additional fields can lead to flattened potentials. 



2 Axions: General Considerations 

The possibility of using an axion in string theory as the quintessence field, with the shift sym- 
metry being broken by non-perturbative effects, has been considered in the past and has some 
well known difficulties. We begin by reviewing this case first. 

Consider a canonically normalized axion field (f) with action, 



S = d^Xy/^g 



(2-1) 



The axion potential in this case is a periodic function with a period of order f a - the axion decay 
constant. For example, a typical form for the potential is 

y = /i 4 cos(£), (2.2) 

where the scale // goes like, 

/i 4 oce- c5inst , (2.3) 

with Si nst being the action of the instanton which gives rise to the potential, and c being a 
constant of order unity. 

For the axion field to play the role of quintessence its potential must be slowly varying on 
cosmological scales and therefore should meet the well-known slow roll conditions. This gives 
rise to the requirements that 

V" r0 _ fv ,s 1 



For a periodic potential, e.g., of form eq. (I2.2[) . it is easy to see that V" JV ~ (V'/V) 2 ~ l//„- 
So the slow-roll conditions give the constraint 

fa > M Pl . (2.5) 
Now axions which arise in string theory typically satisfy the condition, [8|l9|l23~tl2"l] . 

/. < (2-6) 
where is the action of the instanton which gives rise to the axion potential. 
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For the axion potential V in eq.f l2.2p . with eq.f l2.3p . to be of order A 4 , eq. fll.ip . requires an 
action, [8], 

S inst ~ 200 - 300. (2.7) 

Thus we see that string constructions lead to an axion decay constant f a , eq.f l2.6p . which 
is about two orders of magnitude smaller than the value required by the slow-roll conditions 
eq.(l23|). 

With N axions evolving together the mismatch does improve. The slow-roll conditions now 
lead to 

fa > • (2-8) 

However, agreement with eq.f l2.6p . eq.f l2.7p . requires, [H], 

N > 10 4 - 10 5 , (2.9) 

which is a large number indeed. 

In the kind of construction we explore, the key new element is that the resulting axion 
potential is not a periodic function. Instead it is approximately a linear function of the axion 
l 10 l [lj. The potential we get has the form 



v{4>) = A = 7%, (2.io) 

la 



where \l is a mass scale and we have expressed the canonically normalized field <fi in terms of 
the dimensionless axion a as 

= f a a. (2.11) 
The resulting action relevant for the late time evolution of the universe in our model is then 



S axion / d X\/ ~g 



Ml.R-^fUdaf-^a 



(2.12) 



For the potential eq. f)2.10p the slow roll conditions eq. fl2.4p give, 

<p > M Pl . (2.13) 
In terms of the dimensionless axion a this condition is, 

a > %L. (2.14) 

J a 



10 The more correct form of the potential will be described in later sections. The linear approximation is a 
good one when the axion takes values in an appropriate range and otherwise is still a good approximation for 
making parametric estimates, which is our main purpose in this section. 
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The second condition we must impose is that the potential energy in the axion today is of 
order the total energy density in the universe, 



(i a ~ A ~ 10 (eV) (2.15) 
From eq. (12.14ft this gives rise to a condition on the scale of the potential 

^ < T^A 4 . (2.16) 

Mpi 

These are in fact the only two important conditions that need to be imposed. As long as 
a model can be constructed with a linear potential, with a scale fi which meets the condition 
eq. (l2.16p . and in which the axion around the current epoch of the universe meets the condition 
in eq. (l2.14p . one would have a workable model of quintessence. Note in particular, that the 
slow- roll conditions, eq. (l2.14[) . can be met in this model even when f a < Mpi by choosing a 
large enough value of the axion. 

Some additional points are also worth making at this stage. 

First, the rolling axion gives rise to a stress energy of the perfect fluid form with an equation 
of state, 

V = u a p (2.17) 

where 

u^t^^Jpl (2.18) 

2 /2 + V e/3 + 1 K J 



where e is the slow roll parameter, 

Ml, fV'V M% 

pi I — (2.19) 



V 

When <f) ~ O(Mpi) we see that u a can be significantly different from —1, the value for the 
equation of state of the cosmological constant. 

Second, the value obtained for the equation of state parameter (from WMAP+BAO+//o + 
D M +SN) in [25] is 

u BE {z) = w + w 1 ^—, u = -0.93 ±0.13, ui = -0.41+°;£ 2 . (2.20) 
1 — \~ z 

Fitting to the central value in our model leads to the value, 

~ 2.14M Pi . (2.21) 

The time-dependent constraint u>\ is also satisfied. 
Third, the axion satisfies the equation 

4> + 3H(p+V' = 0. (2.22) 
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When the slow roll conditions are satisfied, 

2,H4> + V = 0, (2.23) 

the total change in the axion field during the evolution of the universe, upto the current epoch, 
can be estimated as 

V M 2 

5<j) ~ (PH- 1 ~ —M 2 Pl ~ (2.24) 

v 4> 

For of order the Planck scale we see that the total change in the axion during the evolution 
of the universe is also of order the Planck scale. 

For a workable model, the potential must be slowly varying for this entire range of variation 
of the axion field. In particular, in our case the linear approximation with a slow enough rate 
of variation, should be valid for the entire range of evolution of the axion. This is a significant 
constraint to meet. As was discussed in the introduction it would be difficult to ensure this while 
working within a low-energy effective field theory alone. For once the axion shift symmetry is 
broken one would expect quadratic and higher terms to typically appear 

+ (2-25) 

la m pi 

thereby ruining the flatness conditions required for slow-roll behavior. By embedding this con- 
struction in string theory we can go beyond general field theory considerations though and as we 
discuss below, we will find that in some models these corrections do not arise and the potential 
remains of linear form, with small corrections, for the entire range of O(Mpi) variation in the 
canonically normalized axion. 



3 More Details on the Model 

We now turn to a more detailed description of the model. 

We will work with flux compactifications in IIB string theory. This is a reasonably well 
studied corner of string theory by now, e.g, see [26j|27]. One starts with a Calabi-Yau manifold 
and carries out a suitable orientifolding to allow for the presence of flux. Then turning on flux 
and/or adding branes gives rise to a warped Calabi-Yau internal space. 

The axion field we will consider arises in this setup from the zero-mode of a two-form field. 
There are two-possibilities, the NS two- form, B2, or the RR two-form C2. In the C2 case we 
start with the ten-dim. action, 

.9? 12 

where fi = 0, 1 • • -3 labels non-compact directions and a, b, a'b' label the six compact directions. 
We then reduce to get the four-dim. action. 



(27r) 7 a 



- 2 R- ±-d^C ab d^C a , bl g aa ' g w + 



(3.1) 
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This gives, 

S = d^x\J—g± 



M 2 1 
2 a 2 y ' 



(3.2) 



where the four-dim. Planck scale is, 



2L 6 



M pi = (o M 2 i ( 3 - 3 ) 
and the volume of the internal space is 

V = L\a')\ (3.4) 

with L being the dimensionless modulus. Keeping only the dependence on the overall volume 
modulus and the dilaton in the axion kinetic energy terms we get, 

n 2 M 2 

f2 y S 1VI Pl (o r\ 

Ja JT-- (<3.5) 

As a simple model we can take the internal space to be a six-torus with equal size circles. The 
axion comes from the zero mode of say the T 2 spanned by the first two internal directions. This 
gives, 

a = C 12 (3.6) 

with, 

n 2 M 2 

f 2 _ y S lvl Pl (a 7\ 

In the more general situation of a Calabi-Yau orientifold with several moduli, other moduli 
will also enter in determining f a but the dependence on the overall volume and g s should still 
be parametrically of the form eq. fl3.5p . Since the volume modulus in particular gives some of 
the more stringent constraints we will neglect these additional moduli in our discussion below. 
A more complete analysis will have to include them of course. The case where all the complex 
structure moduli are much heavier and decouple due to a tree-level superpotential, and the 
Calabi-Yau has only one Kahler moduli, will essentially map to the case above *3 

In the discussion above we have considered the case of a RR axion. For an NS axion the only 
difference is that the factor of g 2 on the RHS will be missing in eq. (13.51) . As we will see below, 
it will be advantageous for our purposes to take the axion to arise from the RR field. For this 
reason we mostly present the formulae for the RR case in our discussion. 



11 One situation which can be qualitatively different is if the axion arises from a two-cycle which is localized in 
a highly warped region. In this case f a is suppressed and goes like, f a ~ e Atop , which is roughly the maximum 
value of the warp factor along the localized 2-cycle. As a result meeting the slow-roll conditions requires large 
values of a. We do not pursue this possibility any further below. 
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We see from eg. (13. 5 j) that for g s < L L > 1, / < Mp;, in keeping with the general expecta- 
tions in string theory discussed above o. To meet the slow roll conditions for a linear potential 
we saw above that eq. (j2.14j) needs to be met. From eq. (l3.5p we see that the axion has to satisfy 
the requirement, 

L 2 

a > — . (3.8) 
9s 

It is well known that the NS and RR axions have an approximate shift symmetry in string 
theory. For example the NS two-form has a coupling on the world-sheet, 

< 3 - 9 > 

where A = 0, 1, ... 9 denote all space-time directions. Now for a two- form, Bab, which has no 
field strength, H = dB = 0, it is easy to see that the right hand side in eq. (I3.9j) is a total 
derivative and must vanish in the absence of any boundary for the world sheet. In this way we 
see that there cannot be a potential which arises for an axion coming from the zero-mode of B 2 . 
Similarly, for the RR fields it is well known that the vertex operator corresponds to the field 
strength rather than the gauge potential, thereby again leading to no potential. 

The argument in the previous paragraph is in fact true to all orders in the a' and the g s 
expansion. For the NS case they fail once non-perturbative corrections in a' are introduced 
and world-sheet instantons can give rise to corrections which generate a potential for the axion. 
For the RR case space-time instantons carrying the charges of the relevant Euclidean Dl brane 
are needed. The argument can also fail in the presence of branes. In the presence of a Dp- 
brane for example, it is well known that a constant B 2 field leads to additional charges for 
the brane (corresponding to Dp-branes with smaller p) and correspondingly additional tension. 
This happens in particular for a D5 brane, which in the presence of non-zero B 2 can acquire 
D3-brane charge. By S-duality, and this will be particularly relevant for our model, it follows 
then that this can also happen for an NS 5-brane in the presence of a C 2 field. 

Our basic strategy will be to break the approximate shift symmetry, which prevents a poten- 
tial for the axion, by including branes in a highly warped region of the compactification. This 
will lead to a potential, but one which is suppressed by the warping. 

3.1 The Basic Setup 

Before discussing the breaking of the shift symmetry in more detail let us give some more details 
about the basic setup of the model. This consists of a Calabi-Yau manifold in Type IIB string 

12 More precisely, an instanton contribution which breaks the shift symmetry for a would arise from a EDI brane 
wrapping the T spanning the first two directions. This would have action Si ns t ~ L 2 /g s , so that f a ~ Mpi/S- lns t 
which agrees with eq. (|2.6[) . 

13 More correctly one means B — F where F is the world volume gauge field. 
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theory with three-form and five- form fluxes turned on [26]. An orientifold projection is needed 
to be able to turn on the fluxes, more generally one works with an F-theory compactification. 
The orientifold projection retains states invariant under (— 1) Fl Q<j, where a is a Z2 symmetry 
of the Calabi-Yau manifold. The Kahler moduli arise from non-trivial two-cycles. If there are 
2 non-trivial two-cycles S 1; S 2 then there are correspondingly two Kahler moduli in the parent 
Calabi-Yau. Now if these two-cycles are exchanged by a, then after orientifolding only one 
Kahler modulus, corresponding to the even combination, survives, see [28] for more details. 
Similarly zero modes, related to the 2 non-trivial two-cycles, arise for the B2, C2 two-forms and 
the C4 four- form. For C4, after orientifolding the even combination survives, while for B2, C2 the 
odd combination survives. The orientifold symmetry breaks the supersymmetry down to M = 1. 
The zero mode for C4 and the Kahler modulus which both arise from the even combination of 
the 2 two-cycles give rise to one chiral superfield which we denote by T. The zero modes from 
B2, C2, denoted by b, a, respectively, which arise from the odd combination give rise to another 
chiral superfield denoted by G = b — ra, where r = Cq + i/g s is the dilaton axion field. 

In general there will be several T and G moduli. In the discussion below for clarity, we 
will simplify and only consider the case where there is one pair of two-cycles and therefore one 
resulting T and G moduli each. In this case the real part of T is related to the overall volume 
modulus eq.f l3.4p by 

Re(T) ~ L\ (3.10) 

Besides the various moduli mentioned above additional ones arise from complex structure 
deformations as well. Once the fluxes are turned on these complex structure moduli along with 
the dilaton-axion will acquire a mass at tree- level J26j[29]. In the KKLT model after integrating 
out these moduli there is a resulting superpotential of the form, 

W = W + Ae~ aT (3.11) 

which will stabilize the overall volume and the C4 axion. 

SUSY breaking in the KKLT scenario occurs by adding D3-branes. Let us denote by U mo d 
the scalar potential which is generated for moduli stabilization. Then the SUSY breaking scale 
M$b in the KKLT model is of order, 

M% B ~ U mod . (3.12) 
3.2 Breaking The Shift Symmetry 

We now turn to breaking the axion shift symmetry. In our model the parent Calabi Yau manifold 
has two "throats" , or highly warped regions. Later on we will see that it is safest to have a Z2 
discrete symmetry TZ, present under which the manifold is invariant and which exchanges the 
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two throats. This Z2 symmetry is in addition to the symmetry a involved in the orient if olding, 
and acts non-trivially on the Calabi-Yau orientifold. 

The two-cycles £i,X 2 , descends into each of the two throat. A 5-brane is located at the 
bottom of the first throat, where the warp factor acquires its minimum value, and wraps a com- 
bination of the 2 two-cycles that is invariant under the orientifold symmetry. In the presence of a 
non-zero value for the axion the 5-brane acquires charge, and additional tension, corresponding 
to D3-branes which are filling all the non-compact directions and are point-like in the internal 
space In addition, an anti 5-brane is located at the bottom of the second throat and again wraps 
an orientifolding-invariant combination of the two -cycles. The anti 5-brane can be thought of 
the 5-brane wrapping the two-cycles with opposite orientation. As a result D3-brane charge is 
induced on the anti 5-brane in the presence of the axion and correspondingly it acquires addi- 
tional 3-brane tension. The net configuration is not supersymmetric. For additional details see 
also p. 

Note that the presence of the anti 5-brane is important for charge cancellation. The total 
D3-brane charge must cancel by Gauss's law. Without the anti brane the additional 3-brane 
charge induced on the 5-brane would either not have been allowed at all or at least would not 
be allowed to relax with time. 

For a D5-brane in the presence of B 2 the D3 charge arises due to the WZ term on its world 
volume, 

SwzDb = / B 2 A C 4 (3.13) 



and the additional tension comes from the DBI term, 



S 



DBID5 



^3 / d°xy/-det(g + B 2 ) (3.14) 



(2*) 

For an NS 5-brane, it then follows by S-duality that there is a WZ coupling and DBI term 
involving C 2 

Swznss = / C 2 A C4 (3.15) 



S 



NS5 



^— [ d 6 x^/-det(g + g s C 2 ). (3.16) 
"9s a J 



(2vr) 5 

In the latter case this gives rise to a potential for the axion, 

Here, aa' = j J] C 2 . The factor of 2 is because of including both the 5-brane and anti brane. And 
as in our discussion of the kinetic energy for the axion above, we have only shown the parametric 
dependence on the overall volume modulus and suppressed the dependence on other moduli. 
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The parameter e in eq.f l3.17p arises from the warped geometry. Consider a warped compact- 
ification with metric, 

d s 2 = e 2A{y) dx^dx» + e- 2A{v) g ab dy a dy b (3.18) 
where the warp factor at the location of the 5-brane is e A °. Then e is given by 



e = e 4A °. (3.19) 



It is useful to consider the limit where the axion has a large value, 



9s 

In this case the potential becomes linear 



L 2 

a > — . (3.20) 



V = 2e- — -a. (3.21) 

(2tt)W 2 



Comparing with eq.f l2.10p we see that 

^ = " (O K ,2 ( 3 - 22 ) 

and eq.f l3.19j) implies that 

//oce 4A °. (3.23) 

When eq. fl3.20p is met the 5-brane tension itself is insignificant compared to the contribution 
due to the 3-brane charge. And one can think of the setup as essentially consisting of a stack 
of 3-branes in one throat with another stack of anti 3-branes in the Z 2 image throat. As the 
axion evolves the induces 3-brane charge decreases and along with it the induced tension of the 
3-branes also decreases. 

The linear potential eq.f l3.2ip is a good approximation when a ^> L 2 /g s , but as noted above 
the slow-roll conditions only require that a > L 2 /g s . When this later more general condition 
is met, we must use the full form of the potential in eq. f|3.17|) . Our analysis will be mostly 
parametric in nature and to simplify the algebra we will mostly use the linear form below. 

Two conditions need to be met for a workable model, these are summarized in eq. (12. 14ft and 
eq. (12. 161) . In the limit, eq. fl3.20p we see from eq.f l3.5p that the slow-roll condition eq. fl2.14p is 
met. In addition if the warp factor is small enough the axion energy density 

A 4 (3.24) 



r (2vr)W 2 
can also be of the required small value. 
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We will take the axion a > 0, for simplicity. 
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To investigate this last condition, let us use eq. (13. 3p and express the string scale a' 2 in terms 
of Mpi and the moduli L, g s . This gives 

2L 12 ( A \ 4 

f ~ wy^i [wj (3 ' 25) 

Taking L = 10, g s = 1 and a ~ L 2 and with (A/M P/ ) 4 ~ 10" 123 we get 

e ~ 10" 120 . (3.26) 

While this is a small number the point to remember is that e is determined by the warp factor 
at the bottom of the throat, eq. (13. 19j) and this is often exponentially sensitive to fluxes, as for 
example happens in the Klebanov-Strassler [30] case. Thus modestly small ratios in flux can 
give the required large hierarchy between the energy density in the quintessence field and the 
Planck scale and eq. (13. 24ft can be met while taking the values of L,g s to be quite reasonable. 

In summary, we see that, as a first pass, breaking the shift symmetry by placing 5-branes in 
highly warped throats allows us to meet the requirements for a model of quintessence. 



4 Other Terms in The Potential 

This is by no means the end of the story, however, for we have not included the effects of 
supersymmetry breaking and moduli stabilization as yet. As was emphasized in the introduction 
these are expected to impose stringent constraints on any model. We will now turn to examining 
these constraints within the KKLT model for moduli stabilization 0. The conclusion will be 
that an additional term in the axion potential does indeed arise and typically dominates over 
the one we have kept above. However this term is also linear in the axion and by adjusting the 
warp factor at the bottom of the two throats our model will be viable after all. 

4.1 Other contributions to the Axion Potential 

The first comment to make in studying the interplay of the moduli stabilization potential and 
the axion is that the tree-level potential which arises in the presence of flux in IIB theory does 
not give rise to a potential for the Ci and B2 fields. This is in agreement with our general 
considerations about the shift symmetry not being broken at this order. 

There are two main possibilities to consider next. First, in the KKLT model non-perturbative 
effects are used to stabilize Kahler moduli. These effects might give rise to a potential for the 
axion. Second, such a potential could arise due to the effects of the warp factor. The additional 

15 Another possibility is to use a' corrections for stabilizing moduli [31]. We leave an investigation of our model 
using such a mechanism for moduli stabilization for the future. 
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D3-brane charge that the axion gives rise to in turn back reacts on the geometry and can produce 
a change in the compactification volume. Since the volume modulus has been stabilized such a 
change will come at a cost in energy. We examine these two possibilities in turn below. 

Before doing so, let us ask in general terms when an additional contribution to the axion 
potential, SV can be neglected compared to the term we have already discussed, eq.f l3.17p . 
Clearly the additional contribution should be small compared to the term we keep, 

5V < V. (4.1) 



In addition the added term should not change the validity of the slow-roll conditions, eq.(|27 
If the corrections are polynomial in the axion, this is true automatically once eq. (14. ip is met and 
the slow roll conditions are imposed for the leading term. E.g., one of the slow-roll conditions 
gives 

W\\, 2 [5V\ 2 M 2 Pl ^ 

The slow-roll conditions due to leading term requires that eq.f l2.13p is valid. Eq. f !4.2p then 
follows from eq.f l4.ip . If the corrections are periodic in the axion, the same slow- roll condition 
gives 

f) 2 f«, 

Using eq.f l3.5p we see that this condition is more restrictive than eq.f l4.2p at weak coupling, when 
g s < 1, L > 1. In the discussion below we will find that the corrections which can be neglected 
can be made sufficiently small so that their effects in both eq. (14. ip and on the slow-roll conditions 
will be small. 



4.2 Contributions From Moduli Stabilization 

4.3 Corrections to the Superpotential 

As was mentioned above in a KKLT model type situation the potential for the volume and other 
Kahler moduli arise from non-perturbative effects |12] . Keeping only the overall volume this 
takes the form, 

W np = Ae~ aT . (4.4) 

Corrections to this superpotential which depend on the axion would generate an addition po- 
tential for it. We turn to examining them next. 

The non-perturbative effects responsible for eq.f l4.4p could arise from Euclidean D3-brane 
(ED3) instantons or they could arise due to gaugino condensation in a 3 + 1 dim. non-Abelian 
gauge theory obtained as the low-energy limit on a stack of 7-branes which wrap a 4-cycle in the 
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Calabi-Yau. For the ED3 case it is possible that there are additional contributions which arise 
from bound states of Euclidean Dl branes and D3-branes. These do not seem to be suppressed 
(at large volume) compared to the leading answer above 0. Including them would therefore 
result in a superpotential of the form, 

W = Be~ aT e ClG (4.5) 
where a, c\ are coefficients. The resulting contribution to the axion potential would then be 

SV ~ U mod a ~ M A SB a (4.6) 

which is unacceptably large. 

To avoid this possibility it was suggested in [1] that one consider the other possibility in 
KKLT models and take models where the non-perturbative corrections arise only due to gaugino 
condensation on wrapped 7-branes. This ensures that additional contributions which can depend 
on the axion are highly suppressed. 

The suppression arises due to holomorphy and non-renormalization arguments involving RR 
axions. Let us recount the argument give in [T] here for completeness. The non-perturbative 
superpotential in the Yang-Mills theory is 

W = Ae~ aS (4.7) 

where S is the holomorphic gauge coupling. At large volume 

S = T. (4.8) 

Corrections to eq. fl4.8p which can induce a dependence on the C 2 ,B 2 , axions must vanish at large 
volume, and therefore must be suppressed by inverse powers of Re(T). But by holomorphy this 
means they must also depend on the C4 field which is the partner of the Kahler modulus. Now 
this C4 field is also a type of an axion which arises from the RR sector and such a correction 
term would break the shift symmetry of this field. As we have argued above this cannot happen 
in perturbation theory. This means any correction to eq.f l4.8p must be exponentially suppressed 
in T and this in turn would lead to a further exponential suppression in the superpotential for 
the axion dependent terms compared to the leading contribution in eq. fl4.7p . 

The resulting contribution to the quintessence potential from such a contribution would be 

SV ~ U mod e-? Re ^a. (4.9) 

16 The EDI charge can arise, for example, from induced world-volume flux along a non-trivial two-cycle con- 
tained in the four-cycle wrapped by the ED3. For fixed flux quanta the extra cost in energy for exciting this 
world volume flux vanishes as L — > 00 suggesting that there is no extra suppression. 



17 



Requiring that this correction is subdominant compared to eq. fl3.2ip . gives, 



U mod e-P Re ^a < A 4 (4.10) 

From eq.f l3.12p we get 

e -/me(T) ( _M < 10 -60 (4 n) 

where we have used the fact that M$b > 1 TeV. While the RHS is indeed small the exponential 
sensitivity of the LHS on the volume means this constraint can be easily met for the required 
range of the axion, a ~ L 2 /g s , by taking modestly big values of L and thus Re(T). 



4.4 Corrections in the Kahler potential 

The potential depends on both the Kahler potential and the superpotential so another source 
for additional terms in the axion potential comes from corrections to the Kahler potential. 

In fact the Kahler potential in the tree-level theory itself involves mixing between the T and 
G modulus [28] and as we see below this typically gives rise to an unacceptably big contribution 
to the potential if the axion arises from a B2 zero mode. This is why we are safer in using an 
axion coming from C2, for which no such mixing arises at tree- level, as the quintessence field. 

The tree-level Kahler potential in the case we are considering with one T and G modulus 
is [28], 



K = -3 lot 



T + T+ -e'^Cbb 



(4.12) 



where C is a coefficient determined by the triple intersection numbers of four-cycles. Note that 
only the b field which arise from B2 appears in the Kahler potential which is independent of the 
a component coming from C^- Typically this mixing will give rise to a quadratic term for the b 
axions, 

5V~U mod b 2 . (4.13) 

The slow-roll conditions requires b > L 2 so 5V > M^ B and will therefore be much too big. Since 
the mixing terms do not involve the C2 axions they are safe in this respect. 

The Kahler potential will generically receive corrections beyond tree-level and these can 
break the C2 shift symmetry. Since this symmetry can only be broken through quantum non- 
perturbative effects these effects must involve Euclidean Dl-branes wrapping a two-cycle in the 
appropriate homology class. The resulting correction in the potential will be of the form given in 
eq.( l4.9p and will have a further exponential suppression in the volume. As discussed in eq. fl4.lip 
such corrections can be made small enough by taking a modestly big value of Re(T). 
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4.5 Warping Effects 



Next we turn to warping effects. As we will see these will give rise to significant corrections. 

The essential physics here is that the extra 3-brane charge which arises due to the presence of 
the axion gives rise to additional warping, and this warping in turn changes the overall volume 
of the compactification |32j (for the warping corrections in Kahler potential, see also [331435] ). 
Since the volume has been stabilized already this change comes at a cost in energy which depends 
on the warping and thus the axion. 

A precise calculation capturing this effect is difficult to carry out, at least with our current 
knowledge of flux compactifications. In the KKLT like scenario we are considering here, the 
potential for moduli stabilisation arises from non-perturbative effects. To obtain it one uses as 
input a Kahler potential and a superpotential, with the non-perturbative effects being incor- 
porated in the superpotential. However once the 5 anti 5-brane system are included, SUSY is 
neccessarily broken and it is not so clear if the resulting warping effects caused by the axion can 
be included in this manner anymore. A first principlies calculation of the resulting changes in 
the potential is even more challenging. 

While a precise calculation seems hard, a rough estimate which captures the essentially 
physics above can be made as follows. Let the change in the internal volume, Vj, caused by 
warping due to the axion induced 3-brane charge be 5V}. Then the resulting potential for the 
axion can be estimated to be of order 

8Vt 

5V ~ U mod —^ (4.14) 

VJ 

At first sight one might think that the change in the potential should be quadratic in the change 
in 5Vj and not linear in it, since the moduli have been stabilized and one is expanding around a 
minimum for them. However, a little more thought shows that this is not going to be typically 
true. The reason is that the warping will not affect all terms in the potential in the same manner. 
As a result the potential itself changes once the effects of the warping are included. 



4.5.1 A Subtlety 

Actually there is one subtlety regarding the definition of the internal volume which we should 
address before proceeding. For a warped compactification of the type we are dealing with here 
with metric 

ds 2 = e 2A ^dx^dx» + e- 2A{y) ~g ab dy a dy\ (4.15) 

the internal volume is 




(4.16) 



7 And more generally the size of all four-cycles. 
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However, the four dim. Planck scale, Mp t , when expressed in terms of the string scale and g s is 
not proportional to V/. Rather, due to the warped nature of the geometry, we get 

Mh~4lf d e y^ 4A{v) - (4.17) 
a 9 S J 

This suggests that the natural variable for gravitational purposes, whose fractional change de- 
termines the change in the potential in eq. fl4.14l) is is really Vi defined by 

Vi = J d 6 y^e- 4A{y) (4.18) 

which appears on the RHS in eq.f l4.17j) . Henceforth, this is the assumption we will make in 
computing the corrections to the potential in eq. fl4.14p @. It is worth noting that the warp 
dependent corrections to the Kahler potential calculated for SUSY situations in [32] result in a 
correction are of the type in eq f!4. 14|) . with this definition of Vj. 

Now note that the internal metric in eq.f l4.15p is invariant under a rescaling, e 2A — > X 2 e 2A , 
9ab ^ 2 9ab- While this keeps the volume Vj unchanged it changes Vj. This is because the 
warped metric eq. (I4.15P with this rescaling also changes lengths as measured in the non-compact 
directions, and therefore rescales Mp\. Since we are only interested in using Vj to compute the 
fractional change ^r- in eq. f)4.14p though, any such ambiguity will cancel out in our calculation. 

In fact it will be safest for our purposes to fix this rescaling ambiguity by setting the volume 
of the unwarped metric to be unity in string units, 



d%^g = (a') 3 . (4.19) 

Then working with the resulting expression for e~ 4A , which is now well defined, we can calculate 
Vi and its change. 

To get going first consider a situation where we ignore any effects of warping. e~ iA is then 
a constant, 

e" 4A = C v (4.20) 
The internal volume from eq. (14.161) , eq. (I4.19p , eq. f)4.20p is 

vj = cfV) 3 - (4-2i) 

Using our previous notation, eq. fl3.4p . we get, 

Cf /2 = L 6 . (4.22) 
From the definition eq. (14.181) we then get that 

V I = C l = L\a'f. (4.23) 



18 Taking the fractional change in the volume V[ instead does not change the central conclusions because it 
does not change the dependence on the warp factor e A ° for the resulting corrections. 
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4.5.2 The Leading Effect 



We are now ready to estimate the change in Vj due to warping effects. Consider a simple model 
first of a stack of N D3-branes at some location in the internal space. For simplicity we take 
the internal geometry to be flat, i.e., a torus. The geometry is of the form, 



ds 1 = e^dxpdx? + e~ 2A (dr 2 + r'dQi). (4.24) 
The change in the warp factor produced by the D3-branes is, 

Se-^ = * (4.25) 

where r is a radial coordinate in the internal space measuring distance from the branes and R 
is the AdS$ radius, 

R A = 47rg s Na' 2 . (4.26) 

Since we have set the unwarped volume to be unity in string units, eq. fl4.19p . the radial variable 
r in eq.f l4.24p has a range < r < 0{ya'). The correction to the warp factor are well described 
by eq. fl4.25p for r <C 0{yfa'), when the effects of the "image " stacks needed to implement the 
boundary conditions on the torus are not important. Near the stack of branes, as r — > 0, the 
constant term is relatively unimportant and e~ AA is given by eq. (I4.25j) resulting in AdS^ x S 5 
space. 

The change in the internal volume caused by the warp factor can now be calculated as 

SVj = J d 6 y^6e- 4A (4.27) 

Using, eq. fl4.25p gives, 

SVi ~ R 4 a' ~ 4ng s NL 2 a' 3 . (4.28) 

Now actually in our model what we are interested in is the change in volume caused by the 
axion. The D3 charge that the axion induces is 

AN oc a. (4.29) 

From eq. (I4.23P and eq. (I4.28P we see that the fractional change in Vi caused by the axion is then, 

SVr a 



The resulting contribution to the potential this gives is 



5V ~ U mod g s ^j. (4.31) 
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It is easy to see that this is unacceptably large. Requiring that 



5V < A 4 (4.32) 

gives the condition 

U mod g s ^ < A 4 . (4.33) 
Let us set g s ~ 0(1), with U mo d ~ Mg B , and a ~ L 2 , eq. (I4.33j) then gives, 

L>^~W^V. (4.34) 



A 2 V 1TeV , 

Now L = 10 30 leads to a huge internal space and it is easy to see that the resulting string scale 
is ridiculously low. 



4.5.3 Incorporating The D3-branes 

Our analysis is in fact incomplete for one obvious reason. We have so far only included the 
effect due to the NS5-brane where the axion induces D3-brane charge and tension. In the actual 
situation at hand there is also the anti 5-brane where D3-brane charge and tension is induced. 
Including its effect can cancel the contribution found above to leading order, but not exactly, as 
we will see. 

It is best to work in a situation where the ambient N units of 3-brane charge in the throats 
in which the 5-branes are placed is much bigger than the 3-brane/anti 3-brane charge induced 
by the axion. One can then estimate the effects of the axion to leading order in a/N. The 
presence of the ambient 3-brane charge and related five-form flux breaks the symmetry between 
the backreaction effects of the induced 3-brane charge on the 5-brane and the anti 3-brane charge 
induced on the anti 5-brane, as we will see. 

It is useful to estimate the contributions of the D3-brane charge and tension in two steps. 
The effects of the D3-branes arise from localized sources in the metric and five-form equations. 
The charge of an D3-brane is opposite to that of a D3-brane while its tension is the same. It 
is helpful in making our estimates to think of the AA" D3-branes as a sum of two kinds of 
sources [36J. The first kind is A A" 3-branes with both charge and tension opposite to that of a 
D3-brane. In terms of sources these can be thought of as —AN D3-branes. The second kind 
of source are AA" pairs of D3 and D3-branes, with each pair together having no net charge and 
twice the tension of a D3-brane. We will see that the leading effect arises due to the first kind of 
source and this is exactly canceled, in a symmetric situation, by the contribution coming from 
the D3-branes in the other throat. 

This leaves the contribution from the AA" pairs of 3 and anti 3-branes. The system of D3- 
brane D3-brane pairs placed at the bottom of a KS throat was studied in [T3HT5] . see also [3TH39] . 
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Figure 1: A schematic diagram of the internal space showing the two throats related by the Z 2 1Z symmetry 
containing the NS5-brane and the NS5-brane respectively. 

One important change is that the geometry is no longer of the type eq. (14.15j) with g a b being 
the metric of the Calabi-Yau manifold. Rather the backreaction in this case distorts the metric 
by more than an overall warp factor. The ambient five-form flux, it was found, screens the 
effects of these pairs at large distance. As a result their resulting contribution is suppressed and 
sub dominant. 

Below we first discuss the contribution due to the —AN D3-branes, and then turn to the 
effect of the pairs later. 

Including both the AiV D3-branes located at xo and the —AN D3-branes, which originate 
from the D3-branes and are taken to be at xq, leads to the equation for the warp factor, 

V 2 5e"^ = Cl AN ( S6{X ^1 0) - 5 ^±^l) . (4.35) 

The opposite relative sign mean that they contribute oppositely to the change in the Vi, eq.f l4.18p . 
In general the contributions will not exactly cancel, the residual contribution would then typi- 
cally still be of order eq. (14.31 ft and unacceptably big. However if there is a Z2 symmetry TZ as 
in figure (TJ of the Calabi-Yau space, under which the two throats are exchanged, which is also 
respected by the ambient flux and if this symmetry is then only broken by the brane-anti brane 



pair, then it is easy to see that the two contributions will exactly cance 
for assuming such a symmetry when we described the basic setup r^' 1 31 



This was the reason 
53. 



4.5.4 The Subleading Effect Due to pairs 

The subleading effect is due to the AiV brane-anti brane pairs. The resulting change in the 
geometry was calculated in [13] for UV of the KS throat, and [Hl[l5] for IR. As was mentioned 

19 Ref. also discusses the use of a Z 2 symmetry for similar purposes. 

20 Another possibility, not shown in Figure 1, is that the two throats are themselves contained in a warped 
region - the "parent" throat. In the IR the parent throat splits into the two 5-brane containing throats. 
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it is not simply of the warped form, eq.f l4.15p . Nevertheless an estimate of the resulting change 
in volume can be obtained by just keeping the change in the warp factor. 

In the region far from the tip of the KS geometry where the pairs are located, the geometry 
is essentially AdS§ x S 5 and the perturbation goes like, 

12 

6e~ 4A ~ ANe 4A °R 4 ^- (4.36) 

where the e 4A ° factor arises from the warping at the tip of the throat. Note the l/|r| 8 fall off 
which is much faster than the l/|r| 4 fall-off for a D3-brane source 0. The important point 
here is that the brane-anti brane pairs give rise to a normalizable deformation in the asymptotic 
AdS$ x S 5 region, which falls like l/|r| 8 , this deformation is suppressed by e 4A ° because the warp 
factor at the bottom of the throat suppresses the energy density of the pair in this manner. While 
some of the discussion in [13] is in the context of the KS geometry their main result is more 
general and should apply to other situations as well which are asymptotically AdS§ x S 5 , with 
the warped throat terminating in a region where the warp factor attains a minimum value e A °. 

The resulting change in V/ that eq. fl4.36p leads to can be calculated from eq.f l4.18p . Since our 
purpose is to make a rough estimate we may as well approximate the geometry to be a region of 
AdS$ x S 5 with the AdS$ being cut-off both in the UV and the IR; i.e., of the type considered 
in the RSI model, with metric, 



r 2 D2 

ds 2 = —(dxjxn + — {dr 2 + r 2 dnl) (4.37) 
R 2 r 2 



where the radial coordinate has range r m ; n < r < r max . This metric has the form, eq.f l4.15p . 
with 

tab = $ab (4.38) 

being the flat-space metric. The warp factor at the bottom of the cut-off AdS region is r min /.R 2 
which in this crude model should be equated with e 2A ° at bottom of the KS-like throat giving, 

r 2 min = R 2 e 2A °. (4.39) 



The brane anti brane pair located at the bottom of the throat produces a change in the warp 
factor eq. (14.361) and we interested in the change in V} this leads to. Integrating eq.f l4.18p using 
eq.f l4.38p gives, 

6Vi ~ ANe 4A °R 4 a' 2 -^, (4.40) 

mm 

21 Actually there is a log(r) factor on the RHS but we neglect this below. This factor arises because the total 
D3 charge increases in the KS solution as one goes to large r. For our crude estimate we can neglect this effect. 
Instead we model the total flux in the throat to be roughly constant and of order say N. R which appears on 
the RHS then is related to this flux by eq. (|4.26|) . 
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where we have taken r max — > oo since the integral converges in the UV. Next using eq. (I4.39j) 
gives, 

5Vr ~ ANe 2Ao R 2 a' 2 . (4.41) 

So far we have only estimated the contribution from the region far from the tip of the throat 
where the pairs of 3-branes are located. We can also try to estimate the contribution from the 
region near the tip as follows. Let us continue to crudely model the throat as a cut-off AclSc, x S 5 
geometry, eq. fl4.37p with the pairs located at r min . Very close to the sources, it seems reasonable 
to estimate that the change produced in the warp factor by a D3-D3 pair is of order the change 
produced by a single D3-brane. This results in a change in the warp factor, 

6e~ 4A ~ SN^- (4.42) 

and a change in Vj 

5^ ~ 5Ng s a' 2 R 2 e 2A ° } (4.43) 

where in obtaining eq. (I4.43P this time, since we are interested in the region close to the brane, 
we cut off the integral before r gets too large at r ~ r m ; n . Note, that this result is parametrically 
of the same form as the contribution from the region far from the tip eq. (14.411) . Adding the 
effects of the near-and far regions we would therefore expect an answer of the form eq. (14.411) . 

To complete our analysis we now calculate the resulting term in the axion potential. Since 
AN oc a we get from eq. (l4.4ip . eq. (l4.23p .eq. (l4.14p that the extra potential generated for the 
axion is 

5V~U mod e 2A °^-a. (4.44) 

The estimate for the near-region can be improved in the context of a KS throat using the 
results of [T4"|[T5]. This more careful analysis gives results which essentially agree with eq.f l4.43p . 
eq.( jP3l) . 

In the above discussion we have not kept track of signs and coefficients. The contribution 
due to warping eq. fl4.4ip . we will see below, will typically lead to a contribution in the axion 
potential that dominates over other contributions, including the term due to the induced tension 
of the D3, D3-branes we first considered in eq.f l3.2ip . It will be important for the axion model 
to work that this net contribution is positive. Our analysis above is too preliminary to allow us 
to determine this sign. Physically one would expect the net contribution to be positive, since 
otherwise the added D3-brane charge on the 5-brane anti 5-brane would reduce the energy of 
the system and this energy would decrease as the charge increased. However, one should clearly 
try to determine this from a first principles calculation, along with improving the other aspects 
of this calculation as well. We leave this for future investigation. 
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4.5.5 Final Conclusions 



The changes in Vj, eq.f l4.18p . that arise due to warping could possibly manifest themselves in 
corrections to the superpotential and/or the Kahler potential. For example, Re(T) ~ L 4 in the 
unwarped case and it is more generally related to the size of the non-trivial four-cycle. Thus 
T could change once warping effects are included causing in turn changes in the superpotential 
and the Kahler potential o The warping induced changes might also lead to effects that cannot 
be incorporated into supersymmetric data in such a ready fashion, since the induced D3-brane 
charge breaks SUSY. Either way the resulting corrections to the axion potential will be given 
by eq. fl4.44p . which is the main result from this subsection. 

4.6 An Additional Contribution 

Before concluding this section let us discuss one additional contribution which will turn out to 
be unimportant. 

There is a correction to the axion potential which arises due to the interaction between the 
5-brane and the anti brane in the two different throats. In the limit where the induced D3-brane 
tension dominates over the 5-brane tension this is simple to estimate. The potential reduces to 
that calculated between D3 and D3-branes calculated in [38] . This gives a result, 




where the two throats are located at radial coordinates 7*1,7*2 and for simplicity we have taken 
the internal space to be flat. Here, e = e 4A °, is the warp factor at the bottom of the throat. We 
see that this contribution is suppressed by an extra factor of e compared to the tension term 
calculated in eq. fl3.17p . Since e is very small, as the preliminary analysis which actually under- 
estimates e already found in eq. (13.261) . we see that this contribution will be highly suppressed 
compared to the ones we keep. 




5V 



n - r 2 




(4.45) 



5 The Model: A More Complete Look 



To summaries, we started with a potential of the form, eq.f l3.2ip . 



1 



(5.1) 



V = Cl 



We thanks L. McAllister for a discussion on this issue. 
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Then we found that corrections will generate an additional term, eq.f l4.44p . Using the relation 
eq.f l3.12p this can be written in terms of the SUSY breaking scale M$b a s ; 

V 1 = c 2 Ml B e* A ° ( J^) a. (5.2) 

Here L 6 is the volume in string units, and R is roughly speaking the radius of the AdS-like throat 
regions in which the 5-brane and anti brane's are placed 0. In the supergravity approximation 

R 2 /a > 1. (5.3) 

We also note that e A ° in eq. fl5.ll) . eq. (15.21) . is the warp factor at the bottom of the throat. 
The total potential is then linear in the axion and the sum of the two terms, 

V = V l + V 1 = //a (5.4) 

where the final expression on the RHS can be taken to be the definition of the scale /z 4 . 

We should also note that eq. fl5.ll) and eq. fl5.2p are valid only for a > and also for sufficiently 
large values of a. The more correct form for Vq is in eq. fl3.17l) . we will return to a discussion of 
corrections to the linear form of the potential in §6.3. 



5.1 The Energy Scales 

To understand which of the two terms, eq.( 15.ip or eq. ( 15.21) is bigger let us start by first setting 
g s — L — 1. This also sets 1/a' 2 ~ M Pl . Also we set -j^- = 1. Since we have not been able to 
calculate C2 anyways, we also neglect any dependence on the coefficients ci,C2. Finally we set 
a ~ 0(1) so that a ~ L 2 as is needed for the slow- roll conditions to be met eq.f l3.8p . Setting the 
moduli stabilization contribution V\ to be of order the vacuum energy density, we then get, 

Vi ~ Mg B e 2A ° ~ A 4 . (5.5) 

This gives V to be 

/If 4 A 4 

V ~ ^A<. ,5.6) 

Now if the SUSY breaking scale is Msb ~ 1 TeV, 

M 4 A 4 

Thus Vo ~ V\ and both make roughly an equal contribution to the vacuum energy. However, 
as the moduli stabilization scale and SUSY breaking scale are raised we see that Vq begins 



23 More correctly if the throats are of KS type [30] the five-form flux changes along them and R can be taken 
to be an appropriate average for this radius. 
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to become smaller than V\ and less important. For example, when Msb ~ 10 10 GeV, the 
intermediate scale, Vq is much smaller, 

(y ) 1/4 ~ HT 17 eV. (5.8) 

Let us keep track of the dependence on the volume L and g s also now. We set a ~ L 2 /g s , 
eq.f l2.14p . Setting V\ to be of order the energy density gives, 

Vi ~ Mj B (i?V)e 2 ^ ~ M 4 B (i?V) e 2 ^-^- ~ A 4 . (5.9) 

Solving for e 2A ° gives, 



Using the fact that M Pl ~ L 6 /(ag 2 ) and substituting in eq.f l5.ip then gives, 

^4^'(f)A. (,n, 



L 6 \R 2 J V ^ffs / 

For our approximation of classical supergravity to be valid, g s < 1, L > 1, and i? 2 > a'. This 
means even for C/ mo d ~ (TeV) 4 , Vq < A 4 and thus V will be subdominant. Increasing L,R 2 , or 
decreasing g s will make Vo even smaller. Similarly increasing U mo d will also reduce Vq. 

Our conclusion is that the second term V±, which arises due to the interplay of moduli 
stabilization and warping caused by the axion getting turned on, is typically the dominant 
contribution. We must caution the reader that we have not kept track of numerical coefficients, 
some of these could be large, and our statements here are really only parametric nature. 

The warped down string scale at the bottom of the throats where the 5-branes are present 

is, 

1 An Mpig s a 

m s ~ —=e ° ~ — e °, (5.12) 

\ a' L 6 



where in obtaining the last relation we have used eq. fl3.3p . Supergravity modes (KK modes in 
the throat region) have a mass which is even lower by a factor of \/o^/R where now Rq is the 
radius of curvature at the bottom of the throat F 4 ! . 

For the case we discussed first above, around eq. fl5.6p . with g s ~ L ~ a ~ R \/a' ~ 0(1) 
and Msb — 1 TeV we see that the warped down string modes and KK modes have a mass of 
order 10 _3 eV. As Msb rises the mass becomes even lower. For g s ~ L ~ a ~ R^y/a* ~ 0(1) 
and Msb ~ 10 10 GeV this mass is of order 10 _17 eV. Thus there are many very light particles 
which arise due to string modes and KK modes in the highly warped region. We will have more 
to say about these light particles in §6.1. 



'This might well be smaller than R as defined above. 
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5.2 The Constraints 



It is also worth summarizing the constraints on the various parameters of the model. 

The slow-roll conditions require that the axion satisfy the condition, eq.f l3.8p . Our estimate 
for the V\ term in the potential is valid only when the charge contributed by the axion is a small 
fraction of the total charge supporting each throat, this gives the condition, 

a < N, (5.13) 

Using eq. (13.81) and eq.f l4.26n this give, 

L 2 < R^/a' 2 . (5.14) 

Finally the total volume of the internal space is L 6 (a') 3 , this must be bigger than the volume 
of each warped throat ~ R 6 , leading to 

L 4 > R A /a' 2 . (5.15) 

The summary is that the conditions, 

L 2 < g s a R 4 /a' 2 <C L 4 (5.16) 

must all be met. They are mutually compatible so there is no obstruction to doing so. In 
addition the scale /i 4 defined in eq. (15.41) must meet the condition 

A = A 4 (5.17) 

where A is defined in eq. (11.11) . This does require /i to be a very small energy, but as was 
emphasized in §3.2, where a preliminary estimate for fi was carried out by neglecting Vi, this 
can be arranged by choosing a modestly small ratio of fluxes which results in an exponentially 
small value of e A ° . 



5.3 Additional Comments 

Let us end this section with two comments. 

First, it could be that despite our best attempts at being careful we have missed some impor- 
tant additional contributions to the axion potential. The following general reasoning suggests 
that even if this were true, incorporating such effects would probably lead to a workable model 
within the kind of setup we have explored. One would expect that any additional contribution 
is linear in a/N, as long as a < Af. Also, since the shift symmetry is broken by the 5-branes 
whose effects are suppressed by e A ° , such a term should also be suppressed by a positive power 
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of e A °. If this power is smaller than 2, then this additional contribution could dominate over 
the terms we have kept. However in this case we can simply adjust the warp factor so that 
the resulting scale fi in eq.f l5.4p . after accounting for this term, has the correct value. This will 
typically make e A ° even smaller and thus the particles as the bottom of the warped throats even 
lighter, but, as we will see in §6.1, this is an acceptable price to pay since these particles do not 
lead to any observable signals anyways. 

Second, in the KKLT construction the negative vacuum energy density that arises after 
moduli stabilization is canceled by the addition of D3-branes which are placed in a warped 
throat. It is natural to ask whether in our model one can dispense with these D3-branes and 
their attendant throat, and instead cancel the negative vacuum energy by the axion dependent 
contributions, eq. fl5.4p . 

A simple calculation shows though that this would require too many units of flux N stabi- 
lizing the 5-brane throats. Let the negative energy generated by moduli stabilization be —C, 
including it in the total potential gives, 

V = fi A a-C. (5.18) 
The constant C can be absorbed by shifting the axion 

a^a + Aa (5.19) 

with 

Aa = C/fi\ (5.20) 

The resulting analysis in terms of the shifted axion then exactly agrees with what we had done 
earlier, when C — 0. Since the shifted axion must satisfy the condition, eq. fl3.8p . (to get a model 
of quintessence) equating V with A 4 now will again lead to the conclusion that 

/i 4 < A 4 , (5.21) 

for L > l,g s < 1. The resulting shift is therefore 

Aa = > = 10 60 , (5.22) 

where we have used the fact that C = U mo d > O(TeV) 4 . The linear potential we have used 
though is valid only when eq. fl5.13p is true and this condition involves the total value of the 
axion including its shift. Thus the large value of Aa would require a large amount of flux N, 
which makes this idea implausible. 

Our model therefore must incorporate the two throats with the 5-branes and at least one 
additional throat (or perhaps two which are related by the Z2 symmetry) to contain the D3- 
branes of the KKLT setup. 
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6 Some Features of The Model 



We now turn to discussing some general features of the model. 
6.1 Many Light Particles: 

The model has many light particles. We saw that the warped down string scale at the bottom of 
the two throats where the 5-branes are located is at least 10~ 3 eV and typically is much lighter 
for a scale of SUSY breaking M$b > 1 TeV . Thus there are many very light particles in the 
theory which come from string modes and Kaluza Klein modes localized in the warped region. 

One might be worried that so many light particles would be in conflict with observation. 
However, this is not so. The model is in fact closely related to the Randall-Sundrum II model 
(RS II) [40J in which there is only one brane called the Planck brane. More accurately it is 
akin to a Randall Sundrum I model (RS I) jH] but with the warped down energy scale at the 
"Standard Model" brane being 10 _3 eV or much lower. The RS II model can be thought of 
as a limit of RS I where the standard model brane is moved away to infinity resulting in a 
non-compact dimension. 

We take the standard model degrees of freedom to not live in the two throat regions where 
the 5-branes are located. In terms of the RS model they are then degrees of freedom on the 
Planck brane. The four-dim. graviton is a non-normalizable mode in AdS$ and is localized 
on the Planck brane as well. The first thing to note is that the throat regions make a finite 
contribution, of order R 6 , to the volume of the internal space where R is given in eq.f l4.26p and 
N are the number of five- form flux units supporting the throat. As a result the 4 dim. Newton's 
constant is finite and as we have discussed above can take the required value consistent with 
the other constraints of the model. In addition the couplings of the matter fields with the 
KK modes in the throat are suppressed so that corrections to Newton's gravitational law due 
to KK exchange and corrections to energy loss in gravitational radiation are suppressed by a 
factor of (pR) 2 compared to the leading answers which neglect these KK modes, where p is the 
characteristic energy scale involved. Now in the model the flux units iV and thus R has to be 
big enough to meet the condition eq. (l5.16l) . however this can be typically arranged by taking R 
to be not very different from the Planck scale. As a result this suppression is very significant 
and results in no detectable signal. E.g., taking N ~ 1000, g s ~ 0(1), eq. ( 15.16H can be met 
for L ~ O(10), and a ~ O(100). In this case, R^ 1 ~ 10 16 GeV so that the suppression is of 
order (p/10 16 GeV) 2 and is indeed very significant. Finally as argued in [40] the interactions 
between the four-dim. graviton and KK modes are also highly suppressed, because the four-dim. 
graviton is localized on the Planck brane and a typical KK mode is localized deep inside the 
AdS region. 
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Besides light particles in the two throats, there are also moduli, not localized in the throats, 
which are stabilized by C/ mo( j (e.g., the overall volume). These have a mass 

^moduli ~ -T72-- I 6 - 1 ) 

ivi pi 

For Umod ~ TeV 4 this is also a mass of order 1CT 3 eV. As the scale of moduli stabilization is 
increased the masses of these moduli increase, while, as we have seen above, typically the masses 
of the warped down string states and the KK states move down. 

In constructing a complete cosmological model one would have to deal with the moduli 
problem [12|H3] that arises due to all these light modes. Perhaps after high-scale inflation one 
can arrange so that the KK and string modes in the warped throats are not significantly excited 
and the moduli like the overall volume which are spread out over the whole CY are either heavy 
enough to not cause a problem, or are rendered unproblematic by symmetry considerations, |44j, 
or by thermal inflation |45j. 

In particular it is important to ensure that the two throats are not excited to even small 
finite temperature. This would cause a black brane metric to replace the throat geometry at the 
tip, the two NS5-branes would then fall into the horizon of the black brane and the resulting 
equation of state would be more akin to a thermal gas and unsuitable for quintessence. This is 
especially a concern because of the low warped down energy scale at the bottom of the throat. 

If reheating after inflation results in entropy being dumped into standard model fields, this 
can probably be ensured. The suppressed interactions with the throat excitations will then 
prevent the throat degrees of freedom from coming into equilibrium with the standard model 
degrees as long as the temperature after reheating is somewhat lower than the GUT scale, and 
this can prevent the formation of a black brane horizon 0. 



6.2 Rotation of The CMB Polarization 

It is well known that an axion a, which is a pseudoscalar, can couple to the photon through the 

AL = (aFF (6.2) 

term in the Lagrangian causing the E mode of the CMB to be rotated into the B mode [201146(147] . 
Here ( is a constant coefficient. Such a coupling also causes the direction of linearly polarized 
light to rotate in the course of propagation from a distant source [4"8] . 

The extent of rotation of the E mode to the B mode is parameterized by an angle 

Act ~ (Aa (6.3) 
25 We thank S. Kachru for discussion in this regard. 
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where Aa is the total change in the axion due to its time evolution. 

CMB experiments put a bound on Aa. The bound in [25] is 

Aa = -l.ldeg±1.3deg(stat.) ± 1.5 deg(syst.) (68% CL) (6.4) 

which is currently consistent with Aa vanishing. 

A significant limitation of this paper is that we have not discusses how the standard model 
arises in our construction. As a result we do not know whether a coupling of the type in eq.f l6.2p 
in fact ariseQ- In fact it is more or less clear that such a coupling cannot arise in a SUSY 
preserving manner. For that to happen the SUSY partner of the axion would have to play 
the role of the QED coupling constant. Typically this happens only if the axion partner is a 
non-compact scalar. In our model where the axion arises from Cg, however, its partner arises 
from B 2 and therefore is a compact scalar. 

The coupling to electromagnetism could arise though in a SUSY non-invariant way and this 
could well be allowed if SUSY is broken at a sufficiently high scale. In particular, such a coupling 
would arise if the U(l) gauge group of the SM arose from a D5-brane wrapping the same two 
cycle Hll which gives rise to the zero mode for the axion. On the world-volume of the D5-brane 
would be the couplings, 

1 „ „..„ 1 
2 



S = (2na' 2 )T 5 
Reducing to 4 dimensions gives, 
S = (2vra' 3 )T ; 



-g-F^ + - 



Co A F A F 



(6.5) 



L 



d^x-sf^-F^F^ + -a I FAF 



4 



(6.6) 



Now redefining F^ v so that the gauge kinetic term is canonical we see that the coupling ( in 
eq.flO) is 

C = 7V ( 6 - 7 ) 

As a result 

Aa ~ ^-5a. (6.8) 



Using the observed limits in eq. fl6.4p (converted to radians) we then get 

Aa 9 

T? < 10 



(6.9) 



The canonically normalized field is <fi = f a a ~ a g s Mpi/ L 2 as in eq. (I3.7D . The change in <fi is 
given in eq.f l2.24p . using eq. (I6.9j) this gives, 

M Pl 



< W-% < 10" 



(6.10) 



26 Such a coupling between the axion of interest and SU(3) color gauge field must be absent otherwise the 
resulting QCD induced potential for the axion would be much too big. 
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Or rather the same appropriate orientifold even combination of two-cycles. 
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where we have used the fact that g s < 1. This tells us that to be in agreement with the 
bound eq. fl6.4p the canonically normalized field must have a value which is at least two orders 
of magnitude bigger than the Planck scale. As a result the axion field will evolve very little in 
the course of the universe's evolution. The slow roll parameter e, eq.f l2.19p . for example would 
satisfy the condition 

e < 10- 4 . (6.11) 

This makes the equation of state for quintessence essentially indistinguishable from the cosmo- 
logical constant. 

In summary, we do not know for certain whether the axion couples to electromagnetism in 
the form of eq.f l6.2p . Such a coupling would arise if the U(l) gauge field originates on a D5-brane 
wrapping the same two-cycle from which the axion zero mode also arises. In this case the bound 
on the rotation put by CMB data is very significant. It requires the canonically normalized axion 
field to change very little during the course of the evolution of the universe, thereby making the 
equation of state for quintessence essentially like the cosmo logical constant. In such a situation 
our best hope for telling quintessence apart from the cosmological constant would be to look for 
a signal in the rotation for the CMB polarization itself. 



6.3 The Linear Potential 

The linear nature of the axion potential is valid for an appropriate range of axion values, we 
took this form in our discussions above to simplify the analysis. 

The Vo contribution in eq.f l5.ip comes from the induced 3-brane tension and its more accurate 
form is in eq.f l3.17p with e = e 4A °, as noted in §3.2. As the axion runs towards its minimum at 
a —7- the correction due to the L 4 term within the square root in eq.f l3.17p will become more 
important. 

Also, if we allow for both negative and positive values of the axion, a in eq.f l5.ip . eq.f l5.2p 
should be replaced by \a\. E.g, the more correct form for V\ is 

V x = c 2 Ml B e^[^\a\, (6.12) 

which arises because the change in the volume induced by the warping only depends on |o|. 

With C2 > (which makes physical sense) the total potential Vo+Vi is positive. Now actually 
there is an overall constant C in the potential, related to the cosmological constant, which we 
have not worried about. Including this we get the potential to be, 

Vtotai = V + V 1 + C. (6.13) 

It is worth emphasizing here that the V\ contribution is linear in the axion only in the limit 
when \a\ N where N is the total five -form flux supporting the warped throats where the 5 
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branes are placed. Once \a\ becomes larger there will be corrections and one does not expect 
the potential to remain linear in the axion, for example the warp factor e at the bottom of the 
throat will itself begin to depend on a. As was noted in §4.3 with todays' knowledge of warped 
compactifications we have at best been able to estimate the linear corrections, calculating the 
higher order terms is even harder. For now, we note the possibility that including these terms 
might well give a more rapidly varying axion potential for \a\ > N, and the more rapidly varying 
potential might actually help improve the tracker behavior of this model, to which we turn next. 



6.4 Tracker Behavior 

It is well known that many quintessence models exhibit tracker behavior, see e.g., [3]. The 
tracker solution is an attractor, and at least for some range of initial conditions the system is 
drawn to the tracker solution regardless of initial conditions. 

Here let us consider the linear potential case, with the net potential being given by eq. (15.41) . 
The equation of motion for the axion is given by 

+ 3#0+^ = O (6.14) 

J a 

Let us assume that the universe is radiation dominated, then 

H=i. (6.15) 



The general solution to eq. (I6.14p is 



4 

ci - ^t 2 + czH (6.16) 

5 fa 



where c\, c<i are integration constants. We see that the axion runs to the origin, where its energy 
is minimized in a time determined by c\. 

Let us take the solution with c<i = and some fixed value of C\ and examine perturbations 
around it. There are two kinds of such perturbations. One is time independent and simply 
shifts Ci, the second dies like t~^. Since the first perturbation is constant in time we see that 
this solution is not really a tracker. 

Working instead in an epoch which is matter dominated sets H = 2/3t, the essential features 
of the axion solution found above continue to be the same in this case as well. 



6.5 The Cosmological Constant 

The lack of tracker behavior means that our model does not solve the coincidence problem. The 
initial value for the axion, c\, must be chosen to be just right so energy in the axion field is of 
the right order of magnitude today when the Hubble constant has value Hq. 
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In fact this feature is tied to the issue of the cosmological constant in this model which 
we have mostly ignored so far. Including the constant C in the total potential, eq. fl6.13p . and 
working with the linear potential as an approximation we have, 

Kotai = /i 4 |a|+C (6.17) 

where we now keep track of the fact that the potential really depends on \a\. 

It is now clear that changing the initial conditions for the axion in effect changes C. Another 
way to say this is that shifting a and then changing C appropriately can keep Kotai an d in fact 
the full axion Lagrangian eq. fl2.12j) invariant 0. 

An important idea, for which the landscape of string theory has now provided considerable 
evidence, is that the cosmological constant takes its small value due to anthropic considerations, 
see, e.g., [4"9] - l5"T] . Let us consider how such anthropic considerations might work in our model 
O Let us take the constant C and the initial condition for the axion c\ as the two variable 
which can be varied. One possibility is that the axion sits at its minimum at a = and dark 
energy is entirely due to the cosmological constant with C taking the value A 4 due to anthropic 
considerations. But the more general possibility is that the total energy in dark energy is 
shared more equitably between the cosmological constant and the axion field with both C and 
/x 4 |a| being of order A 4 . In this latter case, the equation of state for dark energy could show 
significant time variation. In fact, knowing nothing better, one would tend to believe that 
this latter possibility is more likely, simply because it is more general. Of course deciding this 
issue more systematically would require a well motivated probability measure in the space of all 
possibilities, which we lack at the moment. 

Let us emphasize that it was important in the discussion of this subsection that the axion 
potential is slowly varying. Instead suppose the scale in the potential eq. (15.41) is, /x 4 ~ M$ B , 
which is the SUSY breaking scale. Then starting with c\ ~ Mpi and taking f a ~ Mpi for 
simplicity, we find that the axion field runs to its minimum in time t ~ Mpi/Mg B which is of 
order the Hubble constant at the time of SUSY breaking. This is too fast to be of any relevance 
today. 



7 Conclusions 

In this paper we have constructed a model for quintessence in string theory. The model is based 
on the idea of axion monodromy. An axion plays the role of the quintessence field, its shift 
symmetry is broken by the presence of 5-branes which are located in highly warped throats. 

28 This is related to some of the discussion in §5.3. 

29 In fact anthropic considerations for quintessence with a linear potential were studied in [17,18 . A probability 
distribution for the parameters of this model, partly based on eternal inflation, was assumed in the analysis. 
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We show that even after including the effects of moduli stabilization and SUSY breaking, the 
resulting potential for the axion can be made slowly enough varying to result in a workable 
model of quintessence. If the canonically normalized axion field has an initial expectation value 
of order the Planck scale, the equation of state of dark energy shows significant time variation 
during the evolution of the universe. 

Our model has many light particles which arise due to the highly warped throats in which 
the 5-branes are placed. The energy scale at the bottom of these throats is at least 10~ 3 eV and 
typically is much lower. The light particles arise from warped down string modes or KK modes 
and are analogous to the light particles in the Randall Sundrum II model which has only the 
Planck brane with a non-compact extra direction. The couplings of these particles to standard 
model fields and to the four dimensional graviton are highly suppressed, making them difficult 
to detect. 

We have not attempt to construct a complete model of cosmology, including for example 
inflation at early times. We have also not attempted to explicitly incorporate the standard 
model fields in the kind of construction we used. This latter issue especially deserves further 
attention because a coupling of the axion to electromagnetism would lead to a rotation of the E 
mode of polarization of the CMB to the B mode, with potential observable consequences. We 
show how such a coupling can arise in our model when the photon is the gauge field living on a 
D5-brane and calculate the resulting rotation effect. 

The cosmological constant question is still left open in our discussion. Anthropic consider- 
ations could well be responsible for its smallness. We suggest that these considerations, in the 
context of our model, would probably lead to the conclusion that the axion field is not exactly 
at its minimum. Rather, it is more likely to be rolling and carrying a reasonable fraction of 
the total dark energy, thereby causing the equation of state for dark energy to vary. Models of 
the type we present here, where the potential for the quintessence field is slowly enough vary- 
ing, might therefore tilt the "betting odds" in favor of quintessence, even though they lack an 
explanation for the cosmological constant along conventional, as opposed to anthropic, lines. 
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